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Abstract 

In this work we develop a general procedure for constructing the 
recursion operators for non-linear integrable equations admitting Lax 
representation. Several new examples are given. In particular we 
find the recursion operators for some KdV-type systems of integrable 
equations. 







1 Introduction 



It is well known that most of the integrable nonlinear partial differential 
equations 



Ut 



admit Lax representation 



F(t, X, It , U x : ' ' ' ; U nx ) ; 



L t = [A,L], 



(1) 



(2) 



so that inverse scattering method is applicable. The generalized symmetries 
M of (ffl) have also Lax representations with the same L-operator 



L tn = [A n , L] , n> 1. 
The recursion operator 71, satisfying the equation (see 

TZ t + [D F ,7Z} = 0, 



(3) 



(4) 



where Dp is the Frechet derivative of the function F, generates symmetries 
of (|IJ) starting from the simplest ones. In general 1Z is a nonlocal operator 
(a pseudo-differential operator). 

The construction of the recursion operator of a given integrable system 
([]]) is not an easy task. Several works are devoted to this subject. Among 
these works most of the authors use @) for the construction of the recur- 
sion operator |§-||. There are several difficulties in this direct approach. 
The main problems are the choices of the order of 1Z and the structure of 
the nonlocal terms. This is an approach, having no relation with the Lax 
representation (Q). 

On the other hand some of the authors used Lax representation for this 
purpose. Most of these works are related to the squared eigenfunctions of the 
Lax operator ||- [IB] and are based on finding eigenvalue equation for the 
squared eigenfunctions of the Lax operator. The operator corresponding to 
this eigenvalue equation turns out to be the adjoint of the recursion operator. 

There is an alternative use of the Lax representation to construct recur- 
sion operators. This approach is based on the explicit construction of the A n - 
operators (^]). It was first used by Symes Adler |Tj| (see also Dorfman- 
Fokas [16 1, Fokas-Gel'fand ]17|]) and Antonowicz-Fordy fl8||,[19|. Although 
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these authors use the Lax representation in different ways their approach is 
basically the same. Symes and Adler use the Gel'fand-Dikii [^(J construction 
of the A n - operators. On the other hand Antonowicz-Fordy determines these 
operators from integrability condition (|3]) and by using an ansatz for A n . 
Their basic aim is to determine the Hamiltonian operators 9\ and 9 2 |21] of 



the equations under consideration. The recursion operator is simply given 
by 1Z = 6 2 0i ■ Their approach is based on some explicit formulas for coeffi- 
cients of y4 n -operator. This is necessary to find the Hamiltonian operators 9\ 
and 02 and it seems that this approach is quite effective to determine the bi- 
Hamiltonian structure for the simple cases but it becomes more complicated 
when L-operator has a sophisticated structure. 

If one is interested only in the determination of the recursion operator 
71, we shall show in this work that it is possible to succeed this without 
any concrete information of the coefficients of A n -operators. We use only an 
ansatz A = V A + R that relates ^-operators for different n. Here V is some 
operator which commutes with the L-operator and R is the remainder. 
We follow this basic idea, partially used by Symes |L|], Adler |fT5|| , Sha- 



bat and Sokolov [] ||22j| , and establish an extremely simple, effective and al- 
gorithmic method for the construction of recursion operators when the Lax 
representation (Q) is given. 

In the next section we consider the case where L is a scalar operator. 
We first consider the case where L is a differential operator and then the 
case where it is a pseudo-differential operator. In each case we present our 
method, discuss the reductions and give examples for illustrations. In section 
3 we consider Lax operator taking values in a Lie algebra. We give our 
method both for the general case and also for the reductions. We give one 
example for each case in the text. Several additional examples are given in 
the Appendices A, B, and C corresponding to all different cases. 

2 Scalar Lax representations. 

First we consider equations with the scalar Lax representations of the form 

1 In 1980 Shabat and Sokolov independently found the recursion operator for the 
Sawada-Kotera equation. This result was published in |22|] . In 22 Sokolov found the 
recursion operator for the Krichever-Novikov equation. 
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L t = [A,L], (5) 

where L is, in general, a pseudo-differential operator of order m and A is a 
differential operator whose coefficients are functions of x and t. 

The different choice of operators A for a given L leads to a hierarchy of 
nonlinear systems (||). It is well known that one can define operators A n by 
the following formula pD[ 



A n = (L"/™) + , (6) 

where L n ^ m is a pseudo-differential series of the form L n l m = Y^oo v ^ anc ^ 
(L n//m ) + = ^2i =0 ViD\ Here are some concrete functions depending on the 
coefficients of L and D is the total derivative with respect to x. 



In p4|,|25l the relationships between the Kac- Moody algebras and special 
types of scalar differential and pseudo-differential operators L were estab- 
lished. All corresponding integrable systems are Hamiltonian ones. For most 
of them a second Hamiltonian structure is not known up to now. 

In this section and Appendices A,B, and C we consider the simplest sys- 
tems from [|4|] and J25| as examples and find their recursion operators. In 
the sequel these examples will be referred as Drinfeld-Sokolov (DS) systems. 
It is interesting to note that in all these examples the order of the recursion 
operator is equal to the Coexter number of the corresponding Kac-Moody 
algebra. 

2.1 Gel'fand-Dikii Systems. 

In this section we shall consider the case where L is a differential operator, 

L = D m + u m „ 2 D m - 2 + --- + u , (7) 

where Ui, i = 0, 1, m — 2 are functions of x, t. In the framework of this 
corresponds to the Kac-Moody algebras of the type A 



m— 1 ' 



To show that (||) is equivalent to a system of (m — 1) evolution equations 
with respect to Ui one can use the following standard reasoning. Set 

L m = (Z/ m ) + + (Z/ m )_, (8) 
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where (L™) + is the differential part of the series L>» and is a series 

of order < —1. Since [L, L™] = we have 

L] = [L, (i™)-]. (9) 

The left hand side of (|9|) is a differential operator, but the right side is a 
series of order < n — 2. Thus both side of (0) are differential operators of 
order < n — 2 and it is equivalent to a system of evolution equations for the 
dependent variables Ui, i — 0, 1, m — 2. This system can be obtained by 
comparing the coefficients of D l where 0, • • • m — 2 in (|3|). 
Since L 11 ^ = LL™ then we have 

A m+n = (LL^) + = L(L^)+ + (L(L^)_) + , (10) 
which leads directly to 

L tn+m = [A n+m , L] = LL tn + [(L (^)_) +! L]. (11) 

The above equation (|ll|) has been given also by Symes fl4| (see also Adler's 
paper ]15[). In his work Symes expressed the coefficients of the both parts of 



(|TTD , in a rather complicated way, in terms of some finite set of coefficients 
of the resolvent for L- operator. That allows him to express L tn+m in terms 
of L tn . This relation gives directly the recursion operator. He gave explicit 
formulas for the cases m = 2 and m = 3. 

In this section we shall show that in order to construct the recursion 
operator it suffices to know only that 

L tn+m = LL tn + [R n ,L]. (12) 
Obviously, it follows from the following 
Proposition 1. For any n 

LA n + R n , (13) 

where R n is a differential operator of order < rn — 1. 
Proof: The relation (|T3|) coincides with (|T0"D if we put 



R n = {L(L&)J) + . (14) 
Since (L™)_ is a series of order < —1, then ord(R n ) <m—l. 
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Remark 1: It follows from the formula 



A, 



(Lm L) + = {L™) + L + L)+, 



(15) 



that 



and 



A n + m A n L -\- R n) 



Lt n L + [L, R n ], 



(16) 



(17) 



where R n is a differential operator of order < m — 1. 

To find the recursion operator we can simply equate the coefficients of 
different powers of D in (0). It is easy to see that in this comparison of 
the coefficients of D l , i — 2m — 2, m — 1 we determine R n in terms of the 
coefficients of operators L and L tn . It is important that the resulting formulas 
turn out to be linear in the coefficients of L tn . The remaining coefficients of 
D l , i = m — 2, • • • , in (O) give us the relation 



( Uq \ 



\ U m- 2/ t 



K 



/ Uq \ 



(18) 



where TZ is a recursion operator. Instead of (12) one can use (0). The 
corresponding recursion operators coincide. 



Example 1. KdV Equation. The KdV equation 



u t = -^(u 3 x + Quu x ), 



(19) 



has representation with 



L = D + u, 



A = (L 3 / 2 ) + . 



(20) 



Since in this case L tn+2 
(II2D takes the form 



= u n+2 and L tn = u tn = u n the main relation 
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u n+2 = (D 2 + u) -u n + [Rn,L], (21) 

with R n = a n D + b n . 

Now if we equate successively to zero the coefficients of D 2 , D and D° in 
above equation, we obtain 

1 3 

a n = -D~ l (u n ), b n = -u n , 

and 

U„+2 = {^D 2 + U+ ^rU x D~ X )u n , 

that gives the standard recursion operator for KdV equation 

U = ^D 2 + u + ~u x D~ x . (22) 

In the same way one can find a recursion operator for the Boussinesq 
equation (see Appendix A). 



2.2 Symmetric and skew-symmetric reductions 
of differential Lax operator. 

The standard reductions of the Gel'fand-Dikii systems are given by the con- 
ditions L* = L or L* = —L. Here * denotes the adjoint operation defined as 
follows. Let L be a differential operator, L = Sa^ D l . Its adjoint L* is given 
by L* = £ (— D) 1 ■ a. t . It is easy to see that if L* = L then m = ord (L) must 
be an even integer. For the case L* = —L it must be an odd integer. 

It is well known that for both reductions all possible A n are defined by (|6|), 
where n takes odd integer values. This condition provides that (A n )* = —A n 
which is necessary for ([jj) to be compatible. 

If L* = L, the formula A n+m = (LL™) + = (L~^~) + gives a correct A n - 
operator since n + m is an odd integer. Thus, in this case Proposition 1 
remains valid and the recursion operator can be found from (12) or (|l 



On the other hand, if L* = —L then both integers m and n are odd and 
hence their sum m + n is an even integer. This means that (L~^~) + can not 
be taken as an A„-operator. In this (skew adjoint) case we must take 
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Ai+2m = {L m )+ = L™) + , 

to find the recursion operator. Following the proof of Proposition 1 we obtain 



Proposition 2. If L* = —L then 

A n +2m = L A n + R n , 



where ord(R n ) < 2 ord(L). It follows from (E3) that 



L 



tn+2n 



L 2 L tn + [R n ,L]. 



Remark 2: Instead of (123D we can use the ansatz 



A 



n+2m 



or 



,4 



n+2m 



L A n L + i? n , 



A n L + R n . 



(23) 



(24) 



(25) 



(26) 



The recursion operators obtained by the utility of (E3|), (|25|), and (p6|) all 
coincide. 

In the works |24j] , |25| more general reductions U = ±L were also consid- 
ered. Here U = KL*K~ Y where K is a given differential operator, such that 
LK~ X is a differential operator. In this general reductions, as well, possible 
^-operators are given by (|6]) with n being an odd integer. The Propositions 
1 and 2 are valid for this general symmetric and skew-symmetric cases and 



hence one can use the equations ([12]), (|2~4]) accordingly to obtain the recursion 
operators. 



Example 2. Kupershmidt equation. This equation 



has the Lax pair 



ut = U5x + I0uu 3x + 25u x u 2x + 20m u x 



L = D 3 + 2uD + u x , A — (Ls' 



(27) 
(28) 



In this case L* = —L therefore we use the equation (|24|) with 



7 



R n = a n D 5 + b n D 4 + c n D 3 + d n D 2 + e n D + f n . (29) 



By equating the coefficients of powers of D in ( 24 ) we obtain 



a 



n 



2 11 1 

~D~ (iin), b n = —u n , c n = -(20uD~ (u n ) + 73u nr 



d n = -(10u x D 1 (u n ) +22uu n + 27u n> 2 X ), 

e n = ^{70u 2x D-\u n ) - 2D-\u 2x u n ) + A0u 2 D-\u n ) - 8D-\u 2 u n ) 
+ 134w n)3:r + 2l2uu n>x + 18Au x u n , 

fn,x = 7^{2Qu Ax D~ l (u n ) + 74u 3x u n + 126u 2x u n ,x + 4:0uu 2x D~ 1 (u n ) 

+ 4Qu 2 x D~ l {u n ) + 13<ou x u H)2x + 27uu x u n + 28u n)5x + <oAuu n ^ x 
+ 16u 2 u n ^ x ), 

and the recursion operator for the Kupershmidt equation: 

U = D 6 + l2uD 4 + 36u x D 3 + (A9u 2x + 36u 2 )D 2 

+ 5(7u 3x + 2Auu x )D + \3u Ax + 82uu 2x + 69t£ + 32m 3 + 

2u x D~ 1 {u 2x + 4m 2 ) + 2(u 5x + 10uu 3x + 25u x u 2x + 20u 2 u x )D~%30) 

2.3 Pseudo-differential Lax operator. 

In this section we generalize our scheme to the case of pseudo-differential Lax 
operators. The only difference is that in formulas like ([H|) and (^) the R n 
operator becomes also a pseudo-differential operator. 

It follows from these formulas that the structure of the nonlocal terms in 
R n is, in general, similar to the nonlocal terms in L since A n j rTn and A n are 
differential operators. 

For skew-symmetric case, A n may be defined by either (|23|) or fl25|) , or 



(|26[). In the pseudo-differential case they are not equivalent in the sense that 
the nonlocal part of R n depends on which ansatz we choose. For illustration, 
let us consider the case L = M D , where M is a differential operator. The 
following lemma shows that if = L or 13 = —L, where 
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l3 = DUD-\ (31) 
then the formulas (p~3|) and ( p5|) are much suitable then (|T6|),(f23|), and (|26|). 
Lemma. Let L' = eL, where e — ±1. Then 

R n = D m - 1 + ■ ■ ■ + oo, /or e = 1, (32) 
where i? n is defined by (0), and 

J R n = J D 2m - 1 + --- + a_ 1 D- 1 , /ore = -l, (33) 
where R n is defined by 

Proof: If L = M D" 1 then L 1 " = e L implies M* = —eM. It is easy to show 
that (L™)' = —L™. Hence (L™)^ = — Ls for an odd integer n. Define now 
a series K n by 

= DK n . 

It is easy to prove that K* = K n . Since K n = (K n ) + + (-Kn)- and (K n )* = 
K n , we have 

(K n y + = (K n ) + , (K n y_ = (K n )_. 

From the last formula it follows that ord(K n )- < —2 that leads to an im- 
portant result 

A n = {L%) + = D{K n )+. 

This implies that 

LA n = M(K n )+, (34) 

is a differential operator. Now using ( |34"D in flT3| ) and ( pop for the cases e = 1 
and £ = — 1 respectively we find the ansatz for A n given by fl32"|) and (53). 



Example 3 (e = — 1). It is known that the KdV equation has, besides the 
standard Lax representation, the following Lax pair: 

L = (D 2 + u)D-\ A = (L 3 ) + . (35) 
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The L-operator satisfies the reduction = —L. According to the formula 
( j33|) we have 

R n = a n D + b n + c n D~ x . 

It follows from (E25h that 



a n = D 1 {u n ), b n = u n , c n = -u n)X -uD l {u n ). 
The remaining equation in fl25|) gives the recursion operator 

TZ = D 2 + Au + 2u x D-\ (36) 



Example 4 (e = 1). DSIII system. The DSIII system ]23J] is given 
by 



"t = -u 3x + 6uu x + 6v x , 
v t = 2v 3x — 6uv x . 

The nonlocal Lax representation for this system is 



(37) 



L = (D 5 -2uD 3 -2D 3 u-2Dw -2wD)D-\ (38) 
A = (L*)+, 

where w; = v — u^x- Since = L we can use ( |3"2D which gives us 

R n = a n D 3 + b n D 2 + c n D + d n . (39) 
By equating the coefficients of the powers of D in (^) we obtain 



dn 



a n = D 1 6 n = 4« n , 

c n = ]^{-QuD- 1 {u n ) + llu n)X + 2D- 1 (uu n ) + 2D- 1 {v n )), 

= ( 6 «2z («n) + 10 Mx Un - 5 « n ,3x + 4 W It^ - 6 V n>x ) 



The recursion operator of the DSIII is found as 
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with 



(40) 



-JZ° Q = D 4 — 8uD 2 - 12u x D - 8u 2x + 16u 2 + 16v + (-2 u 3x + 12 uu x + 

\2v x )D~ l + 4u x D~ 1 u, 
TZl = -10D 2 + 8u + Au x D- 1 , 

Ki = lUv x D + I2v 2x + (Av 3x - UuvjD- 1 + Av^u, 

U\ = -4D 4 + 16uD 2 + 8u x D + 16v + 4u a ,D~ 1 . (41) 

This recursion operator has recently been given in |J . 

3 Matrix L-operator of the first order. 

In this section we demonstrate how our approach, given in the previous 
sections, can be generalized to the case where L is a matrix operator of the 
form 

L = D x + Xa + q{x,t). (42) 

3.1 General Case 

Let us consider the Lax operator (^) , where q and a belong to a Lie algebra 
g and A is the spectral parameter. The constant element a is supposed to be 
such that 

g = Ker (ad a ) © Im (ad a ). (43) 

First let us recall the procedure |24j] of constructing the A- operators for the 
Lax operator (|42|). 

Proposition 3. There exist unique series 

u = -u_i A -1 + ii_ 2 A~ 2 + • • • , Uielm^ada), (44) 

h = h + /i_x A -1 + /i_ 2 A -2 H , hieKer (ad a ), (45) 
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such that 

e°*. (L)=L + [u, L] + i [u, [«, L]} + ■ • ■ = D x + a\ + h. (46) 

Let b be a constant element of </ such that [b, Ker (ad a )] = {0}. It follows 
from (g§) that [6 A n , D x + aX + h] = 0. Hence L] = 0, where 

$ M , = e^" (£>A n ). (47) 
Then the corresponding A-operator of the form 

A M = 6 A™ + a n _! A™- 1 + • • • + a , (48) 
is defined by the formula 

A b , n = ($&,„)+, (49) 

where 

{VL 0O a i \ i ) + = T%a i \ i . (50) 

According to ( f4~7D 

$6,n+l = A $ 6 , n . (51) 

Hence 

= (A $ 6 , n )+ = A ($&,„) + + (A ($6,n)-) + - (52) 

The last formula shows that 

Afe.n+l = A4, n + -Rn, -Rn £ 9, (53) 

where i? n does not depend on A. Substituting (|53|) into the Lax equation 
L tn+1 = [A bin+1 ,L] we get 

L tn+1 =\L tn + [R n ,L}. (54) 

Using the ansatz (|54D one can easily find the corresponding recursion opera- 
tor. 

Example 5. The system 
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1 2 

U t = -- U xx + U V, 

is equivalent to the nonlinear Schrodinger equation, has a Lax operator 

The Lie algebra g in this example coincides with sl{2). 
Using fl5"4|) with R n = ( fln ^ n J we find that 

a n = - (vu n + uv n ), 

h - I -_i 



and the recursion operator of the system (pq) is given by 

_ f—\D + uD ~ l v uD~ x u \ , , 

3.2 Reductions in matrix case 

In the general case considered in the previous section the A n -operators belong 
to Lie algebra 

a+ = {Sfc=o «i A*, Oi e <?, «G Z + }, (58) 
that is a sub-algebra of the Lie algebra 

a= {T, K _ rx> a i X i , caeg, keZ}. (59) 

A standard cr-reduction is defined by any automorphism a of the Lie 
algebra g of finite order k. Because a K = Id, the eigenvalues of a are e l ,i — 
0, • • • ,k — 1 where £ is primitive K-root of unity. 
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Let gi be a eigenspace corresponding to eigenvalue e % . Then the following 
reduction a,j G ft, where i = j (mod k) in fl58|) and ( |59|) is compatible with the 
equations ([3]). Note that according to this definition a G ft and the potential 
q(x, t) in (|42|) belongs to ft or , the same, satisfies <r(g) = q. 

It is easy to see that, to satisfy such a reduction, we must use the ansatz 

Ab tn+K = X K A b:n + R n , (60) 

where 

Rn = r K _ 1 A K " 1 + --- + r , r< G ft. (61) 

Further generalizations are associated with modifications of sign " + " in 
( pDD which corresponds to the simplest decomposition of algebra a into direct 
sum of two sub-algebras 

o = a+©a_, (62) 

where a + is given by (|58|) and 

a_ = {EZ^ a i A i , G </}. (63) 

The sign " + " in (|50D is the projection of onto a+ parallel to a_. If we have 
a different decomposition ( |62]) then the construction from Proposition 3 is 
also valid but we have the following condition 

R n Ea + n\a-, (64) 

instead of R n G g. If we also have the cr-reduction we must use the most 
general ansatz (^) where 

i? n Ga + nA K o_. (65) 

Example 6. Let us consider the following equation 

1 3 3 3 

u t = - u xxx u xx u H — u u xx uu x u, (66) 

4 8 8 8 

where u is a square matrix of arbitrary size, or more generally, u belongs to 
arbitrary associative algebra /C. This equation has a Lax representation with 
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Here 1 is the unity of /C. The reduction ( p7|) can be described as follows 



see [26|). The Lie algebra 5 is the algebra of all 2 x 2 matrices with entries 



belonging to JC. The automorphism a is defined by 

a(X)=TXT-\ (68) 

where T = f J ^1 ) Obviously a 2 = Id and eigenvalues of a are 1 and 
— 1. The corresponding eigenspaces are 



and therefore the coefficients a; in (|59|) have the following structure 



a 2i = ( n ° J » a 2i+i = ( ° * ) ■ (70) 



*y ZJ ~^ v* 

The sub-algebra a + is given by ([58]), where the coefficients have the structure 
( |70|) and, additionally ao = f g |j ^ . The sub-algebra a_ has the following 
form 



tL. = X3?^a i A 4 , (71) 



a 

where ao is of the form ao = [ ] , a £ K, 



a 

The ^-operator for (|66"D is given by formula A = ($^ 3 ) + (see ( ^9|) ), where 
J J j and " + " means the projection onto a + parallel to o_. 



According to (65), R n is of the form 



It follows from 

L tn+2 = A 2 L t „ + (73) 
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that 



&n,x [Oji) ^] &n C-n 0, C n b n Cln,x 0, 
dn b nx ^ &n 0, ~t~ Cn,x Cn H 0, 

U n +2 — —dn,x + [^nj U] . 

Finding a n ,b n ,c n , and d n from this system we obtain the following recursion 
operator 

U = —(D + ad u ) (-D + R u ) (2D + ad u y l (D + L u ) D (2 D + ad,,) -1 , (74) 

where i? u and L u are the operators of right and left multiplications by u, 
respectively. 

Note that in the commutative case (|66|) coincides with the modified KdV 
equation. It is easy to verify that (|74|) becomes the standard recursion oper- 
ator of modified KdV equation. All factors in (|74]) have to be regarded as op- 
erators acting on (non-commutative) polynomial depending onu,^, u xx , 

4 Conclusion 

In this work we devoted ourselves in the construction of recursion operators 
when the Lax representation is given. We have shown that our approach 
can be easily generalized to all cases where L-operator is a polynomial of A. 
It would be interesting to generalize it for the cases of more complicated A 
dependence of L as well as for the cases 2 + 1 dimensional equations, Toda- 
type lattices and ordinary differential equations. 
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5 APPENDIX A 

In this Appendix we give an additional example to Section 2.1. 



A. Boussinesq System. 

The Boussinesq equation 

u u = -^{uax + 2(u 2 ) 2x ) (75) 
can be expressed in the form of a pair of first-order evolution equations 



u t = v x , 

v t = --(u 3x + 8uu x ). (76) 
This system has a Lax pair 

L = D 3 + 2uD + u x + v, A=(Li) + . (77) 

To construct the recursion operator for this system we use the equation (|12]) 
with the differential operator 

Rn = a n D 2 + b n D + c n . 
By equating the coefficients of the powers of D in we find 

2-1 1 
a n = -D (u n ), b n = -(5u n + D (v n )), 

I 

C n = -(6f„ + 8MD _1 (M n ) + 10u ntX ), 

y 



and after that we obtain the recursion operator of the form ( 10 ) for ([73) with 

TZ° = 3v + 2v x D~\ 
Tl\ = D 2 + 2u + u x D-\ 

1 A a 10 o 16 9 ,2 16 

= -(-D i + —uD 2 + 5u x D + 3u 2x + —u 2 + (-u 3x + —uu x )D x ), 

ft* = 3v + v x D-\ (78) 
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6 APPENDIX B 

In this Appendix we give additional examples to the Section 2.2. 



B.l. Sawada-Kotera Equation. 



The Lax pair for the Sawada-Kotera equation [27 



u t = u 5x + 5uu 3x + hu x u 2x + 5u 2 u x , (79) 

is given by 

L = D 3 + uD, A=(L?)+. (80) 

In this example, = —L, where V 1 = D~ X L* D and L is skew- symmetric, 
then we use (p4|). The operator R n has the same form as ( p9| ) with the 
coefficients given by 



1 5 1 

a n = -D~ l (u n ), b n = -u n , c n = -(5uD~ l (u n ) + 29u niX ), 

d n = - (5u x D _1 (w n ) + lAuu n + 26u n , 2x ) , 
y 

e n = 27 ( 10w 2x ^"V™) - ( M 2x«n) - D~ l (u 2 u n ) + 5M 2 D _1 (M r , 
+ 28u n! 3 X + 32uu n , x - 32u x ii n ) , f n = 0. 

The recursion operator is given as 



U = D 6 + 6uD 4 + 9u x D 3 + (9u 2 + llu 2x )D 2 + (10u 3x + 2luu x )D 
+5u ix + 16uu 2x + 6-u 2 + 4-u 3 + (u 5x + 5uu 3x + 5u x u 2x + 

5t?w :r )Zr 1 + w :c Zr 1 (w 2 + 2w 22; ). (81) 

B.2. DSI System. 

The DSI system 
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u t = 3vv x , 

v t = 2v 3x + 2uv x + vu x , (82) 
has a Lax representation with 

L = [D 3 + (u + v) D + X - {u + v) x ] [D 3 + {u-v)D + ^{u- v) x ], (83) 

A = (d) + . 

Here R n is a differential operator of order 5 and since L is symmetric we again 
use the equation ( |T2"D . The expressions for the coefficients of the operator 
R n are very long and complicated. Hence we do not display them here. 
We find that the recursion operator TZ of this system is of the form (|40|) 
where 



U° = -AD 6 -2AuD i -72u x D 3 + 2(-A9u 2x - 18u 2 + A2v 2 )D 2 

+ 10(-7u 3x - 12uu x + 30vv x )D 

- 26u 4x - 82uu 2x - Q9u 2 x + 222vv x + Ulv 2 x - 16m 3 + 4Sv 2 u 

+ 2(— 2u 5x — 10uu 3x — 2hu x u 2x — 10u 2 u x + 15v 2 u x + 30vv 3x 

+ 45v x v 2x + 30wu x )D -1 + 2u x D~ 1 (3v 2 - 2u 2 - u 2x ), 

Tl\ = 168t»D 4 + 20Av x D 3 + 6(2lv 2x + 32uv)D 2 + 6{A0vu x + 7v 3x + 22uv x )D 

+ 6(13f u 2x + lQu x v x + v 4x + 5uv 2x + Avu 2 + 12v 3 ) 
l^vv^^v + 2u x D- 1 (Quv + 9v 2x ) 



Kl = 56vD 4 + 268v x D 3 + 2(2A3v 2x + 32uv)D 2 

+ 2{3Qvu x + 2\9v 3x + lQQuv x )D + 2{27vu 2x + 92u x v x + 99v 4x + 99uv 2x 

+ Avu 2 + 12v 3 ) + 2(10vu 3x + 35u 2x v x + A5u x v 2x + 10uvu x + I8v$ x 

+ 30uv 3x + 10u 2 v x + ^v^D" 1 + 2v x D' 1 (3v 2 - 2u 2 - u 2x ), 

K\ = 108D 6 + 216uD 4 + A32u x D 3 + 6(81u 2x + 18u 2 + 22v 2 )D 2 

+ 6(A5u 3x + 36uu x + 70vv x )D 

+ 3(18u 4x + 18uu 2x + 9u 2 x + 98vv 2x + Q7v 2 x + 32uv 2 ) 

+ 3Q(2v 3x + 2v x u + vu x )D~ l v + 2v x D- 1 (Quv + 9v 2x ). (84) 
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B.3. DSII System. 

The DSII system fl§, [[25| 

u t = 3v x , 

v t = -2(v 3x + uv x + vu x ), (85) 
has a Lax representation with 

L = {D 6 + uD 3 + D 3 u + {v + ^u 2 )D + D{v + ^u 2 ), (86) 
A = 

Since L is symmetric we again use the equation (|T2"D . In this case the operator 
R n is given as follows 

R n = a n D 5 + b n D 4 + c n D 3 + d n D 2 + e n D, (87) 

where 

1 1 5 

c„, = l[5uD- 1 (u n ) + 3D- l (v n ) + 29u n , x ], 
y 

d n = ^ [5w x D" 1 («„) + 26u n , 2x + I4:uu n + I2v n ] , 

e n = ^[5(2u 2x + u 2 + 3v)D' 1 (u n ) - 3D~ x (vu n + uv n ) + 9uD -1 (f„) 

- 2D' 1 (u 2x u n + ^m 2 m„) + 54u x u n + 28u Ht3x + 32(uu n>x - u n u x ) + 42tvJ . 
The recursion operator (f40|) for the system can be found as 

-R° = -D 6 -6uD 4 -9u x D 3 ~ (llu 2x + 9u 2 + 42v)D 2 

+ (—10m 33; — 2\uu x — 84w x )D — 5u ix — 16uu 2x — 6u 2 — 60v 2x — 4m 3 — 2Avu 
+ ( — — 5uu 3x — 5u x u 2x — 5u 2 u x — 15vu x — 15v 3x — Ibuv^D^ 1 
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- u x D 1 (2u 2x + u 2 + 3v), 

11° = -42L> 4 - A8uD 2 - 87u x D - 6(7u 2x + u 2 - 6u) + 27v x D~ x - 3u x D~ x u 

n\ = 28vD A + 106v x D 3 + (165v 2x + 32uv)D 2 

+ (54vu x + 132v 3x + 74v x u)D 

+ 30vu 2x + 79u x v x + 54t) 4x + 57uv 2x + 4u 2 v — 24v 2 

+ (10vu 3x + 25v x u 2x + 30u x v 2x + 10uvu x + 9v$ x + I5uv 3x 

+ 5u 2 v x — 15ff x )D _1 — v x D (3v + u 2 + 2u 2x ), 

K\ = 27D 6 + 5AuD 4 + 135u x D 3 + 3(5Au 2x + 9u 2 - 22v)D 2 

+ 3(36u 3x + 27uu x — 28v x )D + 3{9u^ x + 9uu 2x + 9u 2 x — 21v 2x — I6vu) 

- 18(v 3x + u x v + v x u)D~ 1 -3v x D~ 1 u. ( 

B.4. DSIV System. 

The DSIV System p^j , |25|j which is also known as the Hirota-Satsuma 
system 

1 

u t = 2 U3x ^ UUx ~ ® VV X7 

v t = -v 3x ~ 3uv x , (89) 
has Lax representation with 

L = (D 2 + u + v)(D 2 + u- v) , A=(lA) + . (90) 

Since the operator L is symmetric we use the equation fll2|). I n this case the 
operator R n has the same form as ([39]) with coefficients given by 



1 7 1 

a n = ~D~ l (u n ) , b n = ~u n --v n , 

c n = \[QuD~ l {u n ) + 2D~ 1 {uu n ) - AD~ l {vv n ) + YIu n , x - V2v n , x \ 
8 

d nx = -}-[6u 2x D~ l {u n ) - l2v 2x D~ l {u n ) + 30u x u n - 8u x v n 
16 

+ 2Auu n , x + 15u n>3x - 12v x v n - 8uv ntX - 20vv niX - 28v nj3x \. 
The recursion operator (f40|) for the given system is 
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K°q = ^D 4 + 2uD 2 + 3u x D + 2u 2x + A(u 2 - v 2 ) 

+ (3uu x - 6vv x + ^u^D' 1 + u x D~ l u, 

H\ = -5vD 2 - Av x D - v 2x - 4:iiv - 2u x D~ 1 v, 
5 

%l = ~2 V x D ~ 3 y 2x - ( v 3x + 3uv x )D~ + v x D^ l u, 

TZ\ = —D 4 - AuD 2 - 2u x D - Av 2 - 2v x D~ 1 v. (91) 

B.5. N=3 Hirota-Satsuma System. 

This system is given by [Q 

u t = -u 3x + 3uu x + 3(-v 2 + w) x , 
1 

v t = -^v 3x -3uv x , 

w t = -~ w 3x - 3uw x . (92) 

This is an example for N = 3 system which covers some other N = 2 systems 
as special cases. For instance letting w = we get DSIV and v = we get 
DSIII systems. 

The corresponding Lax pair is 

L = (D 2 + 2u — 2d) (D 2 + 2u + 2v) + Aw, A=(L%)+. (93) 

In this case the operator L is symmetric and hence R n has the same form as 
391) with the coefficients 



7 

a n = D -1 («n), b n = -u n + v n , 

c n = ■^\2uD~ 1 (u^) + AD~ l (uu n + w n - 2vv n ) + I7u n , x + 12v n>x ], 
d n ,x = ^[l2u 2x D~ 1 (u n ) + 2Av 2x D~ 1 (u n ) + 60u x u n + 16u x v n + I5u nt3x 

o 

+ A8uu n>x + 2Av x u n - A0v x v n + 20v n>3x + 16vv n;X + 20w n>x \. 
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The recursion operator is given by 



(K K K \ 
n= \ n l n\ n\ , (94) 

where 



■Rl = -D 4 + AuD 2 + 6u x D + A{u 2x + Au 2 - Av 2 + Aw) 

+ A(-u 3x + 3uu x — 6vv x + 3w x )D~ x + Au x D~ x u, 

H\ = -2(5vD 2 + Av x D + v 2x + 8uv + Au x D~ l v) } 

Ul = 5D 2 + 8u + Au x D-\ 

nl = -5v x D-6v 2x -2(v 3x + Qv x u)D- 1 + Av x D- l u, 

n\ = -D 4 -8uD 2 - Au x D + 8{8w - 2v 2 ) - 8v x D~ l v -8D~ x w x , 

nl = 4(^zr 1 + 2ir 1 u c ), 

nl = — 5w x D — 6w 2x — 2(v 3x + Qw x u)D~ l + Aw x D~ l u, 

nl = -16vD- 1 w x -8w x D~ 1 v, 

n\ = —D 4 - 8uD 2 - Au x D + 16(w - v 2 ) + Aw x D~ l + \QvD~ 1 v x . (95) 

APPENDIX C 

In this Appendix we give an additional example to Section 3. 

C.l. Non-Abelian Schrodinger equation. 

This is the system given by 



u t 



-- u xx + UVU, 



■ V TT — VUV, 



(96) 



where u and v belong to K, (see Example 6 for the notations)). The Lax 
operator of ( |96| ) is given by 
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The corresponding formula (|54D reduces to 

C +1 y)^c («o 



where 



i? n = _^ V (99) 

The formula (p^) gives us both a n ,b n , c n and the recursion operator 7?.. They 
are given by 

«n = + "n), K = ^U n , C n = ~U n , (100) 



^ = 1 f-D + R u D- l R v + L u D- 1 L v R u D' 1 L u + L u D~ l R u 
2 y — L v D 1 R v — R v D 1 L v D — R v D 1 R u — L v D 1 L u 

(101) 

C.2. Non-Abelian Modified KdV Equation. 

The standard non-Abelian Modified KdV equation is given by 

13 3 
u t = -u xxx - -u x u 2 - -u 2 u x . (102) 

The Lax representation of this equation is given 

The recursion operator TZ can be found from (|60| ) and (plf). In our case the 
automorphism a is the same as in Example 6 and formulas (|60"D and (|6l|) 
give us 

° M «+^-A 2 (° U A+[R n ,L], (104) 



v n+ i J \v n 
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where 

*=tt 0)^(0 ly <«*> 

Using ( |104j ) we find a n ,b n ,c n , d n from the following: 

bn ^ n ^Ti) ^rijX ^ W Q. n -|- C n d n 0, 

bn,x ~t~ U ~\~ U b n -\- d n C n 0, d n x -\- C n ^x \pn ^] 3 

= d>n,x + Ki 11 ]' 

The resulting recursion operator is given by 

TZ=^(D-ad u - D- 1 ■ ad u ) (D - {L u + R u ) D~ x (L u + R u )). (106) 
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